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In this paper we suggest a simple mathematical procedure to derive the classical 
^ ■ probability density of quantum systems via Bohr's correspondence principle. Using 

Tij- . Fourier expansions for the classical and quantum distributions, we assume that the 

(N 

Fourier coefficients coincide for the case of large quantum numbers n. We illustrate 
the procedure by analyzing the classical limit for the quantum harmonic oscillator, 
although the method is quite general. We find, in an analytical fashion, the classical 
distribution arising from the quantum one as the zeroth order term in an expansion in 

X 

\ powers of Planck's constant. We interpret the correction terms as residual quantum 

effects at the microscopic-macroscopic boundary. 
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In physics, a new theory should not only describe phenomena unexplained by the old 
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theory but must also be consistent with it in the appropriate hmit 111 . In this sense, newto- 
nian mechanics can be recovered from relativistic mechanics in the domain of low velocities 
compared with the speed of light in the vacuum. Since its formulation, quantum mechanics 
has established itself as the most successful physical theory for the description of microscopic 
systems, such as atoms and elementary particles. Unlike special and general relativity, rela- 
tions between classical and quantum mechanics are more subtle, given that the conceptual 
framework of these theories are fundamentally different. While in classical mechanics it is 
possible to know the exact position and momentum of a particle at any given time, quantum 
mechanics only specifies the probability of finding a particle at a certain position j2|. 

The first statement of a mathematical procedure to obtain the classical limit of quantum 
mechanics can be traced back to Max Planck 3|]. He postulated that classical results can 
be recovered from quantum ones when Plancks constant is taken to zero. Planck originally 
formulated this limit to show that his energy density for black body radiation approaches 
the classical Rayleigh- Jeans energy density when ^ — t- 0. A different approach is due to 
Niels Bohr ^. He postulated that the classical behavior of periodic quantum systems can 
be determined when the principal quantum number is large. Bohr enunciated it in this way 
because in his model of the hydrogen atom the transition frequency between two neighboring 
energy levels tends to the classical orbital frequency of the electron when n ^ 1. Some 
researchers, however, have argued that the two methods are not equivalent 5|-[3]. 

Textbooks and articles on quantum mechanics usually discuss a variety of ways to make 
the connection between classical and quantum physics. Most of them are based on either 



Plancks limit or Bohrs correspondence principle. For example, the WKB SHlQl] and quantum 
potential [ll| methods and the phase space formulation of quantum mechanics are discussed 

nn 

using Planck's limit, while some authors |2|, 112| compare the classical and quantum proba- 
bility densities for both position and momentum, showing that these distributions approach 
each other in a locally averaged sense (coarse-graining) for large quantum numbers n. There 
are other proposals, like Ehrenfests Theorem 13|] , based on semi-classical approximations to 
quantum mechanics. Another method is by means of coherent states. The standard coherent 



states of the one-dimensional harmonic oscillator |14|-|l6l| are localized wave packets which 
follow the classical equations of motion. However, for non-quadratic Hamiltonians this only 
holds approximately over short times. 

Wigners phase-space formulation of quantum mechanics offers a comprehensive frame- 
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work in which quantum phenomena can be described using classical language. The Wigner 
distribution function (WDF), however, does not satisfy the conventional properties of a 
probability distribution [l7|; e.g., WDF is in general positive semi-definite. Therefore, in 
order to interpret it as a classical probability distribution, strictly one needs to restrict the 
analysis to situations where it is non-negative (this is the case for coherent and squeezed 



vacuum states only) 18|, ll9|. W. B. Case has made a careful discussion of the classical limit 



and its difficulties via WDF 



20|. 



According to Bohrs correspondence principle, classical mechanics is expected to be valid 



in the regime in which dynamical variables are 



units 211 . In addition, some authors 



m 



22 



arge compared to the relevant quantum 
23[ suggest that we must compare the same 



physical quantities in both approaches , e.g. probability distributions and not trajectories 
or wave functions. 

In 1924, Heisenberg made an attempt to give Bohr's correspondence principle an exact 
mathematical form in order to apply to simple quantum systems. He suggested that for 
a classical quantity f{t) in the case of large quantum numbers, the following approximate 
relation is valid: 



{i^n+m \ f \ i^n) = {n + m\ f \ n) exp[i/h{En+m - En)t] ^ exp[imu{n)t] 



where fm{n) is the mth Fourier component of the classical variable / and u{n) is the 



classical frequency 2J, |25| . The application of this procedure, however, was limited to the 



study of light polarization in atoms subject to resonant fluorescence [26|, 

In 1926, E. Schrodinger proposed a different application of the correspondence principle 
applied to the quantum harmonic oscillator. His approximation consists of adding all the 
wave function oscillation modes, generating a semiclassical wave packet 28|], from which 



other interesting ideas have recently evolved 29|, |30| . On the other hand, dis crep ancies and 



discussion remains about the adequacy of Bohr's correspondence principle |3lN34j. Some 
suggest that the harmonic oscillator does not have a true classical limit when described by 
means of stationary states 35| and others argue that this system violates Bohr's correspon- 
dence principle [36 1. 

In this paper, we suggest a conceptually simple mathematical procedure to connect the 
classical and quantum probability densities using Bohr's correspondence principle. 
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It is well know that for periodic systems, the quantum probability distribution (QPD) 
pQ'^i (^x,n) is an oscillatory function, while the classical probability distribution (CPD) 
p^^{x) does not have this behavior. However, both functions can be written as a Fourier 
expansion, i.e. 



p'^A^ (x, n) = j /«^^ (p, n) e^'-^dp (2) 
P^^(x) = J fHp)e''^dp (3) 

where (p, n) and Z*"^ (p) are the quantum and classical Fourier coefficients, respec- 
tively. In addition, we know that for simple periodic systems these distributions approach 
each other in a locally averaged sense for large quantum numbers. This implies that the 
Fourier expansion coefficients should approach each other for 1: 



f^"{p.n)^f^^{p) (4) 

the principal quantum 



Note that the 



In order to make this comparison we first substitute the value of 
number n by equating the quantum and classical expressions {2], 12 
Planck constant keeps a finite value, so ^-dependent corrections may arise, as implied by 
equation(jl]). 

Our proposal can be summarized as follows. First we calculate the coefficients of the 
expansion (p, n) by using the Fourier transform of QPD, and then obtain its asymptotic 
behavior for large n. We then equate the classical and quantum expressions for the energy, 
to define the value of the principal quantum number. Finally calculating the inverse Fourier 
transform we obtain, at least in a first approximation, the CPD. The procedure can be also 
applied to probability distributions in momentum space. 



We now apply this procedure to the harmonic oscillator 
dimensional harmonic oscillator is given by 



21 



22]. The QPD for a one- 



p^*'(a;,n) 



(5) 



where a = One of the main differences between the classical and quantum descrip- 
tions of the harmonic oscillator is that the QPD is distributed completely throughout the 
X-axis, while the CPD is bounded by the classical amplitude. However, when increase the 
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value of the principal quantum number n, the QPD exhibits a confinement effect, akin to 
the classical behavior. 



We now calculate the Fourier coefficients. 



many handbooks of mathematical functions 37 



le corresponding integral can be found in 



38|: 



\2muhJ 

where L„ is a Laguerre polynomial of degree n. We remark that the mathematical 
structure of the coefficients f^^^ (p, n) is similar to the Wigner function for the harmonic 



oscillator 



39| . but formally different, due to the dependence of the wave functions on parity 



[40|. Technically, the WDF is a member of the Cohen class of phase-space distributions which 
is related to the fractional Fourier transform 4l| , and not with the usual Fourier transform 



as is the case for the expansion coefficients. 

The asymptotic behavior of Fourier coefficients for n large is also well known. Szego 
finds the following iterative relation: 



42 



71 

2 



F (u^) = e-'^Ln {u^) ~ Jo (2VNu 

t^F {t^) Jo (2VNuj Yo (2VNt^ - Jq (2VNt^ Yq (2VNu 



where Jq and Yq are the usual Bessel functions of the first and second kind respectively, 
and = n+^. Szego shows that in — oo limit the iteration terms are strongly suppressed 
compared to Jq Bessel function. 

2 

Using the above relation and = , we can write the asymptotic expression for the 
Fourier coefficients as follows 



TT 

2 



X 



y/2mujh 



t^F {t^) 



Jo(^]Yo(2VNt 



Jo(2^^t)YJ^^ 



dt (8) 



\ n / \ / \ J \ n / - 
Finally, we compute the inverse Fourier transform. The first term can be obtained di- 
rectly, while the iterated terms can be written as dimensionless integrals 



1 1 y./ Tl\kfh\ 



2k 



(9) 
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where S = nmuxQ is the classical action and the ik = {x,Xq) is the kth dimensionless 
integral. In particular: 

H (x, xo) = / tiae*°^ / /3Vo [Jo (a) ^ (/?) - Jo (/?) ^ («)] (10) 







We can also evaluate higher order iterations in a simple fashion 42|. 



Note that the first term in Eq. ([9]) is /i-independient and corresponds exactly with the 



CPD 



Il2| . The remaining terms are proportional to increasing powers of 4, which are very 



s 

small for classical systems, so these terms are strongly suppressed compared with the CPD. 
A residual oscillatory behavior, as observed in the QPD is preserved through the harmonic 
behavior of the iterated integrals. If we now consider Plancks limit, the classical result is 
exactly recovered. This, however, is not necessary, as the correction terms are very small 
and seem to reflect a residual quantum behavior at the classical level. 

A complete agreement of both the position and momentum distribution functions at the 
classical limit is necessary for the theory to recover the classical results in the appropriate 



energy limit [43|. In this case, due to the symmetry of the harmonic oscillator, the QPD in 
momentum space can be obtained easily, so the asymptotic behavior of the QPD for large 
quantum numbers is given by: 

2 

where po is its maximum momentum, S = is the classical action and ik{p,Po) is the 
same dimensionless integral defined by Eq. ( |T0|) . 

Expectation values of physical quantities can be calculated using our previous results and 
the classical values are then recovered, i. e. 



~ 4l (12a) 
{f)-vlL (12b) 
h) ~ EcL (12c) 



where we have not included the correction terms. These results do not ensure that the 
time dependence of position and momentum operators defined by the Heisenberg equation 
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reduces to the classical equations of motion, due to the fact that the classical limit is not a 
single trajectory, but an ensemble of trajectories. 

To summarize, the classical limit problem has been debated since the birth of quantum 
theory and is still a subject of research. In this paper, we present a simple mathematical 
formulation of Bohr's correspondence principle. We consider the simplest quantum system, 
the harmonic oscillator, and obtain exact classical results. We believe that this approach 
illustrates in a clear fashion the difference between Plancks limit and Bohrs correspondence 
principle. 

Finally, using this simple procedure we find corrections to the exact classical result as 
a series in the ratio ^, which is very small for classical energies but not zero. It would be 
interesting to test whether this energy dependence could be observed for the case of real 
quantum systems approaching the microscopic-macroscopic boundary. 

We thank Alejandro Frank Hoeflich and Jose Adrian Carbajal Dommguez for his valuable 
contribution to the fulfillment of this work. 
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